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Abstract 

In this article we consider one dimensional model of an ultra short 
pulse propagation in isotropic dispersionless media taking into account 
a nonlinearity of the third order. We introduce a method for Maxwell's 
equations transformation based on a complete set of projecting operators. 
The operators generally correspond wave dispersion branches. As a sim- 
plest result of the method application we derive a system of equations 
describing dynamics of ultrashort pulses of opposite directions of propa- 
gation. We show that in such way the generalized Short Pulse Equations 
of Shafer and Wayne is obtained if the only directed wave is initialized. 
The effects of the unidirectional pulses interaction are traced. 

1 Introduction. 

Recently the trend towards shorter pulses duration account was successfully 
realized in a waverange 680nm to infrared [T]. The resulting equation (named 
as Short Pulse Equation - SPE) was properly investigated numerically 2 and 
it was proven to be integrable O |4] . The derivation of the SPE was based on 
few approximations, included the dispersion relation adjusted to the mentioned 
waverange, third order nonlinearity account and unidirectionality of a pulse 
propagation realized via slow time variable introduction in the frame following 
the pulse. 

A significant interest to the both directions of propagation account arises 
in different context, e.g. in resonator theory. One of examples of such account 
based on the second order wave equation is given in Ref. [5] . Author shows that 
unidirectional approach for hybrid electromagnetic fields allows to write down 
first order wave equation for pulse propagation. Moreover he shows in his article 
a generalization of undirectional ^ and bidirectional [7] approach presented in 
previous works. 

Based on the interest and its wide development we suggest a general method 
[SJin] implementation adjusted to the problems to be considered and which links 
unidirectional pulse propagation approach j^Oj with derivation of SPE \X, and 
its generalizations. 

In Sec. II we outline an idea and realization of projection operators that 
specify evolution operator subspaces. Later on we present an application of this 
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approach to the derivation of bidirectional wave system for hybrid fields H] 
and in final section our results show in what conditions a transition to equations 
of earlier works [I] are realized. 



2 Basic theory 
2.1 The model outline 

Our starting point is the Maxwell equations for non-magnetic medium in the 
Lorentz-Heaviside's unit system 

V-D = 0, (1) 
V-B = 0, (2) 

VxE . --f. 

c at 

VXH . 122 ,4) 

c at 

completed by material relations within the choice 

D = E + 47rP, H = B. (5) 
The polarization vector P have the form 



'NL- 



(6) 



We restrict ourselves by one-dimensional model as one of Ref. [T] , x-axis chosen 
as a direction of a pulse propagation. We assume that = and = Q 
taking into account the only polarization of electromagnetic wave. This allows 
us to rewrite the Maxwell equations as 



1 dDy _ dB, 



c dt dx 
I dB, _ dEy 

c dt dx 



where Dy ^ Ey + Airx'^^^Ey + (PNL)y, 



In such model, implying the properties of isotropic opticaly non-active media 
[TTl [12], the third order nonlinear part of polarization supposed to have the 
form: 

(P^L),=4.(xf4^,^). (8) 

These results will certainly simplify the further derivation procedure in the part 
of nonlinearity account. However up to this point we are unable to tell whether 
the unidirectional approximation would become more or less robust for descrip- 
tion of pulse propagation. To establish a point of reference we provide steps 
which lead to SPE equation through application of projection operators. 
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2.2 On projection method. 

As a first step to illustrate the method [5] we define projection operators for a 
simple case of linear isotropic dielectric media where Dy = eEy and e — l+Awx^- 
In this case we can rewrite Eq. ([7]) as matrix equation of a type 



tpt = Lip, 



where the field vector 



enter the matrix operator equation as 



(By), )-[cd._. )[By )-[ c{E^), ) ■ ^^"^ 

Consider now a Cauchy problem for the system Eq. (jlOp . Applying Fourier 
transformation on x and hence rewriting this equation in frequency domain, one 
arrives at ordinary equations system 

{ioByik,Lu)) J \ kc{E,{k,uj)) )■ ^ ' 

Let us next search for such matrix Pi that Pi 5* — '^i and P2^ = ^2 to be eigen- 
vectors of the evolution matrix in Eq. Moreover, the standard properties 
of orthogonal projecting operators 

P,*P, =0, Pf = P,, = 1 (12) 

are implied. For Eq. ffTTI) Pi get the form 

p^ = \i \ ^ V P2 = \i V V (13) 



2VV^ 1/ 2V-V^ 1 , 

Performing the inverse Fourier transform yields the x-representation of the op- 
erators [Sj.In this simplest case we consider the matrix elements of projecting 
operators do not depend on fc, hence its x-representation coinside with the It- 
representation Eq. 

With the projection operators we can introduce new variables A and 11 



which correspond to left and right direction of wave propagation [S]. Moreover 
comparing new variables A and 11 to variables presented by Kinsler et. al. [TU] 
ours have similar form to their ones. The differences are caused by the fact that 
we have taken into account linear interaction between pulse and media. However 
our form of Eq. (jl4[) is exactly determined by dispersion relation uj = ^ from 
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Eq. PT|) . What is more this allows us to present in simple way both electric 
and magnetic field. [3] 

E,=A + n By = V^{A-n). (15) 

This correspondence Eq. (jl4ll5p is one-to-one local map and hence allows to 
determine initial conditions in the Cauchy problem for both left and right wave 
variables (A,n). It also gives a possibility to follow waves, extracting data 
in each time t. The example we consider is simple, but contains all principle 
ingredients of the method [8" . More complicated example of the next section 
show what changes if the projection operator matrix elements depend on k. 



2.3 Dispersion account: unidirectional waves subspaces. 

Within this formalism we reproduce calculations presented in ^J going down 
to a unidirectional waves subspace. If next one assumes that medium of prop- 
agation can be modeled as made of free atoms interacting with the external 
electromagnetic field, one can present an expression on x^^^ of the form |12j 



X 



(1) 



n- \- (16) 



Then for silica fibers and for light in the visible to mid-infrared range there are 
three resonances of importance which occur at wavelengths of A = 0, 068 . . . fim, 
A = 0, 116 . . . fim and A = 9, 896 . . . fim. With those values we have approxi- 
mated x*^^-' by 

(1) _ _(1696A2_ 0.4079A^ 0.8974A^ 
^ ~ A2 - (0.0684)2 + A2 - (0.116)2 ^ A2 - (0.986)2' ( > 

where A is the wavelength expressed in fim. Following [T] we focus on the 
propagation of light in the infrared range with A = 1600 — SOOOnm. In this 
range we can approximate 

X^'H^)^xi'^+xi'^X'- (18) 

If we plug Eq. ^TE\i in e and to Eq. (fT5|) then the projection operators get the 
form 

H Vl + 4-X^^^+4-X^'^# 1 /' 

with the dependence on k that arises from deispersion account. Similar form 
has the second operator 

To get the x-representation for the result we expand the coefficients in Taylor se- 
ries and perform inverse Fourier transformation. With new projection operators 
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we rewrite relations Eq. (|14p in new conditions 



A = VI + 4^X^ + ^^^==(*ax)"^i?, + I , (21) 
^ 2^/1 + 4^x1,'^ 



V 2^/1 + 4^x1,'^ 

3 Nonlinear ity account: general dynamics equa- 
tions, SPE revisited 

3.1 Nonlinear terms as perturbation 

If we account the nonlinearity Eq. ([8]) then Eq. (ITOl) get the form 

d 

— *-L* = N(*). (22) 
Simphfying Eq. dH) to 47rf ^x^^^^^^ [IlIIS] 



d f Ey\_( ld,\fEy\ . d f ix(3)ij3 



dt\B,_ ) \cd, J \ J- ^""dt \' )■ ^^^^ 

Applying Pi Eq. (fT3| on the LHS of Eq. ((23)) and respect the projectors 
property [P^, - L)] = 0, yields 

|--l)p,* = P,N(*). (24) 



dt 

Finally, present the result 



Repeating our calculations from Eq. ((22)) to Eq. ((25)) with use of second projec- 
tor operator Eq. ((20)) than we will obtain a system of equations which describes 
interaction between two waves propagating in opposite directions. The system 
of equations have the form 

\ n£Vin\ ix(^)(A + n)3 



iVin ) - I -c^n j = 2^^^ I -^x(^)(A + n)3 ^ 

^A \ / f^v/iA\ , / ix(^)(A + n)3 \ ' ^ ^ 
cV j-2-^Ux(^)(A + n)3 ) 



where e = 1 + 47rxo^'' + 47rx2^'c^(j9a;)~^ is the operator in x-representation, 
which inverse and other functions we understand via corresponding Taylor series 
expansions. The couples of equations in Eq. (|26p are in fact equivalent. 
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3.2 The generalized Shafer- Wayne SPE 

Let us choose the first equation from the second system. 

= 2^^x(3)(^^n)3. (27) 

If we mean a pulse launched from a right end of a fiber we can consider the only 
direction of propagation. Hence we have chosen A as a dominating left wave. 
Assume therefore that 11 = and differentiate twice with respect to x produces 
the generalized SP equation 



Following [1] we can make an multiple scales ansatz 

A{x,t) = {KAoi4',xi,x2) + K^Ai{4>,xi,x2) + ...) (29) 

where 4> = and Xn — /^""x. If we consider all terms 0{k'^) then we reach the 
final result 



- ( 3cV 1 + 4.xr - 2(1 + 4.x^-) ) - I 4.xr - ^^^^^^ 

= 2-x(^^^(Ao)^ 

At this point we obtain a second order differential equation which describes 
ultra short pulse launched in one direction. To fulfil this description we have 
to present initial-boundary conditions, that is immediate corollary from the 
projection fixed by Eq. (I?!]) within the choice 11 = 0. It reads as the correlated 
action of magnetic and electric fields. 



4 Discussion and conclusions 

As we derive the system of equations Eq. (1261) describing interaction between 
two waves propagating in two directions, we would like to remark that such 
interaction is weak in case of a long optical fiber exited from both ends. However 
the significance of this phenomena is obvious: it give a new possibility investigate 
a nonlinearity and measure nonlinear constants. The fundamental importance 
of the opposite waves interaction phenomena account is quite clear in a case of 
optical resonators, when the effects of interaction are cumulated. Going down 
to unidirectional case we acquire a SPE equation with modified coefficients 
within the demonstration of the general method fSl. One of the purposes of 
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this publication is to attract attention of researches in optics and other fields 
to the projecting operators method. Some interesting features of the projection 
operator method applications to optical problems can be also seen in the paper 
of Kolesik et. al. jl4j. In acoustics there are many interesting applications 
and important development with nonlinearity account in projecting operator 
recently published in Ref. |15) with references therein. 

We have also demonstrated how to obtain SPE with a projection operators 
method for unidirectional pulse propagation in isotropic, dispersionless media 
containing nonlinearity of the third order. The achieved results gives generalized 
SPE and correction to coefficients in it. 
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